
capture log close
log using $data/EDS_11_mobility_patterns, text replace

clear
clear matrix
clear mata
set mem 2g
set matsize 800
set more off
set maxvar 20000
set trace off
set varabbrev off


#d;
	  global cov "_Isu_schule_1 _Isu_schule_2 _Isu_schule_3 
	_Isu_ausbila1 _Isu_ausbila2
	_Isu_kat_al_2 _Isu_kat_al_3 _Isu_kat_al_4 
	_Ichild_1 _Ichild_2
	_Isu_eintri_1 _Isu_eintri_2 _Isu_eintri_3 _Isu_eintri_4 _Isu_eintri_5
	_Isu_eintri_6 _Isu_eintri_7 _Isu_eintri_8 _Isu_eintri_9 _Isu_eintri_10 _Isu_eintri_11
	_Ibefstat_1 _Ibefstat_2 _Ibefstat_3 _Ibefstat_4
	_Idauloint__8 _Idauloint__9 _Idauloint__10 _Idauloint__11 _Idauloint__12 
	_Idauloint__13 _Idauloint__14
	female east_germany
    german migration_background 
    married  
	benefit_status benefit_level log_benefit_level
    lifetime_ue_age lifetime_empl_age last_wage last_wage_mis 
    region_ba_aloq_zivall gdp_real_pc region_ba_vacancies 
	sector_manufacturing sector_service  log_mobirate_avg5
	openness_st conscient_st extraversion_st neuroticism_st loc_st
	write_ger speak_ger write_eng speak_eng
	homeowner  
	   ";		

	   
		#d;
		global outc"
	    w2_relocation50
		w2_commute50
		";
	  #d cr   
	   

		use $data/data_pair, clear


		
		bys pers_id: gen insample = 1 if _n==1
		
		
	
		
		
		
	  foreach type in  $outc {
	 	
			
		
		qui reg `type'  log_intensity_alt  $cov county_pair_*  [pw=1/fweight_reg], cluster(region)
	  
	  matrix b=e(b)'
		matrix v==e(V)
		matrix se=vecdiag(cholesky(diag(vecdiag(v[1..1,1..1]))))'
		count if e(sample)==1 & insample==1
		matrix n=r(N)
		sum `type' if e(sample)==1 & insample==1
		scalar m=r(mean)
		scalar r = 100*(b[1,1]/m)
		matrix d_`type'=b[1..1,1..1],se,r
		
		
		
		
	
}/*type*/
#delimit;
	matrix d=
	 d_w2_relocation50\
	 d_w2_commute50
		   

	;
	
	matrix rownames d = 
	"Relocation 50km+" 
	"Commuting 50km+" 
	;
	#delimit cr

	matrix colname d = "Coef" "SE" "%-effect" 
	mat list d 
	estout matrix(d, fmt(4)) using $data/Table_A10_survey.tex, style(tex) replace
	
	
	
log close	

